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We present an approach to spectropolarimetry which requires neither moving 
parts nor time dependent modulation, and which offers the prospect of 
achieving high sensitivity. The technique applies equally well, in principle, in 
the optical, UV or IR. The concept, which is one of those generically known 
as channeled polarimetry, is to encode the polarization information at each 
wavelength along the spatial dimension of a 2D data array using static, robust 
optical components. A single two-dimensional data frame contains the full 
polarization information and can be configured to measure either two or all of 
the Stokes polarization parameters. By acquiring full polarimetric information 
in a single observation, we simplify polarimetry of transient sources and 
in situations where the instrument and target are in relative motion. The 
robustness and simplicity of the approach, coupled to its potential for high 
sensitivity, and applicability over a wide wavelength range, is likely to prove 
useful for applications in challenging environments such as space. ©2012 
Optical Society of America 

OCIS codes: 120.5410, 260.5430. 



1. Introduction 

The polarization of light provides a versatile suite of remote sensing diagnostics. In astron- 
omy, polarization is used to study the Sun and Solar System, stars, dust, supernova remnants, 
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and high-energy extragalactic astrophysics [1]. The astrophysical mechanisms by which po- 
larized light is produced range from scattering phenomena to the interaction between high 
energy charged particles, and magnetized plasmas. Beyond astronomy, polarization is used 
in remote sensing, medical diagnostics, defense, biophysics, microscopy, and fundamental 
experimental physics, e.g. [2]. 

Accurate, precision polarimetric methods usually require rapidly modulating, often fragile, 
parts and are inherently monochromatic, e.g. photoelastic modulators (PEMs), ferroelectric 
liquid crystals or liquid crystal variable retarders (LCVRs) in tandem with phase locked 
photomultipliers, or synchronized charge shuffling on a charge-coupled device (CCD) detector 
for area detection [1]. Lower accuracy techniques typically require sequential measurements 
of the target using rotating waveplates and polarization analyzers. Here we describe a method 
to encode polarimetric information over a wide spectrum in a single data frame, using static 
optics. This approach alleviates errors introduced by the need to match sequentially acquired 
data, and eliminates the need for fragile or rapid modulation, yet may be able to accomplish 
high accuracy, precision measurements. The methods, of course, have their own implicit 
sensitivities and concerns, as we discuss below. 

A particular interest of the authors, which serves as a useful illustrative example, is the 
use of precision circular polarization spectroscopy as a remote sensing biosignature and a 
potentially valuable tool in searches for biological processes elsewhere in the Universe. The 
circular polarization spectrum is sensitive to the presence of molecular homochirality, a strong 
biosignature, through the combined optical activity and homochirality of biological molecules 
[3,4]. Biologically-induced degrees of circular polarization have been found in the range 10 -2 
to 10~ 4 for a variety of photosynthetic samples, with an important correlation between the 
intensity spectrum and polarization spectrum [3]. Hence, precision full Stokes polarimetry 
and wide spectral coverage are required. Furthermore, the target scene and instrumentation 
may be in rapid relative motion, compounding the difficulties of acquiring the data using 
traditional polarimetric techniques. A large number of photons must be accumulated in a 
short period of time. The techniques presented in this paper may provide a means to make 
this type of polarization measurement, in addition to providing a robust method for acquiring 
less precise spectropolarimetry in a straightforward fashion. Furthermore, the approach is 
applicable across a wide wavelength range, and, as well as in the visible, can work equally 
well in the ultraviolet, where for example chiral electronic signatures are generally strongest, 
to the infrared, where polarimetry goes hand in hand with probes into the geometry and 
physical characteristics of dusty regions of the universe. 

A variety of similar concepts are available under the generic title of "channeled polarime- 
try" [2]. These typically fall into two classes: channeled imaging polarimetry (CIP) and 
channeled spectropolarimetry (CS), following terminology of [2]. To simplify, the CS meth- 
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ods typically encode the polarization information as an amplitude modulation directly on 
the spectrum, derived from a polarization optic whose retardance is a function of wavelength. 
As an example, the spectral modulation principle for linear spectropolarimetry [5] can reach 
a precision of at least 2 x 10 -4 [6]. The CIP methods, by contrast, use a polarization optic 
whose retardance is spatially varying, so that the polarization information is encoded as 
a set of spatial fringes onto an image [7]. These two approaches, as well as a number of 
technical issues that arise in each case, are described in some detail in [2]. Previous authors 
have used multiorder retarders, birefringent wedges, pairs of birefringent wedges, and Savart 
plates individually or in combination for these two applications [5,7-16]. Typically, the po- 
larization information is extracted from the data using Fourier methods. Another approach 
to single-shot imaging polarimetry and spectropolarimetry is the wedged double Wollaston 
device, which yields multiple images on a detector with polarization axes at different angles 
and allows retrieval of the Stokes parameters through combinations of the images [17-19]. 

The approach explored in the current paper is to disperse the spectral and polarimetric 
information along two orthogonal directions, a "spectral" dimension for the spectroscopy and 
a "spatial" dimension for the polarimetry. The amplitude modulation of the encoding of the 
polarization information is independent of the choice of spectral resolution. The two aspects 
of the measurement, the spectroscopy and the polarimetry, may be optimized independently. 
The complete spectropolarimetric information is encoded on a single data frame, and may 
be derived using straightforward analytical techniques. 

Poisson photon counting statistics play a critical role in astronomical polarimetry. To 
measure a polarization degree of 10 _ri , it is necessary to collect (at least) 10 2n photons. 
For example, to measure p ~ 1CT 4 , it is necessary to accumulate 10 s photons. A typical 
astronomical CCD has a well-depth ~ 10 5 electrons per pixel, requiring 10 3 pixel readouts. If 
the data are needed in, say, 1 s in one pixel, this multi-readout approach becomes prohibitive. 
A solution is to spread the illumination across many pixels, as is done for high signal-to- 
noise-ratio photometry with the Hubble Space Telescope [20]. Making a virtue of necessity, 
if we use optics which spread the light of a spectrum perpendicular to the spectrum, then 
we can exploit the width of the broadened spectrum to encode the polarimetric data. 

Sections 2 and 4 discuss a variety of configurations that accomplish this goal. Sec. 2 
starts with linear polarization (equivalently any two of the Stokes polarization parameters), 
followed by a discussion of our analysis methods in Sec. 3. Configurations that enable full 
Stokes spectropolarimetry are presented in Sec. 4. Sec. 5 describes practical implementation, 
sensitivities, and an approach based on calibration. Sec. 6 provides an example application. 
Finally, we make some conclusions in Sec. 7. 

The different embodiments of the underlying approach described in Sees. 2 and 4 highlight 
different aspects of the method. In the end, we anticipate that the most useful realizations of 
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the concept will be the double wedge for linear polarimetry, Sec. 2B, and the double-double 
wedge for full Stokes spectropolarimetry, Sec. 4.B.3. The other subsections introduce new 
ideas incrementally, while these two sections capture the final products for the two types of 
polarimetry. 

We use the conventional Stokes vector formalism to quantify the polarization of light with 
S = (I, Q, U, V) where / is the total intensity; Q, U decribe the linear polarization and V 
the circular polarization. The normalized Stokes parameters (q, u, v) = (Q, U,V)/I represent 
the fractional polarization state. The degree of polarization is given by p = ^ q 2 + u 2 + v 2 
and the direction of linear polarization by ip = | tan -1 (w/g). 

2. Concept for Linear Polarimetry 

We envisage a spectrum of light broadened in a direction orthogonal to the dispersion di- 
rection and sensed using a two dimensional area array such as a CCD. This broadening can 
be spread along a segment of a conventional long slit spectrograph, for example, with the 
length of the entrance slit providing the spatial dimension in the detected two dimensional 
spectrum. To introduce amplitude modulation along the slit (x direction), we introduce a 
retardance gradient 4>(x) along x using a birefringent wedge (or wedges) followed by a polar- 
ization analyzer, such as a dichroic polarizer or polarizing prism (see Figs. 1 and 2). It would 
be possible to carry out the polarization analysis immediately in front of the detector array, 
as in [7]. However, performing the polarization analysis as early as possible in the optical 
path yields better robustness against polarization introduced by the instrumentation optics. 
Furthermore, the polarization optics can be more compact and easier to characterize, since 
light from all wavelengths is analyzed using the same optical elements. Hence, we prefer to 
insert the polarization optics immediately adjacent to the spectrograph's entrance slit, and 
allow a long-slit spatial segment to project through the instrument to become the detector 
array's spatial dimension. 

2. A. Single birefringent wedge 

To lay the groundwork, we initially consider just a single birefringent wedge. The wedge 
thickness gradient is oriented along the slit, while its fast axis is oriented 45° with respect to 
the slit. The analyzer's transmission axis is parallel to the slit, though it could alternatively 
be orthogonal to it. If we define the Stokes Q direction as also being parallel to the slit, then 
we can consider a uniformly illuminated slit with the beam entering the slit, orthogonal to the 
slit plane. If the incoming light is polarized with its electric field along the slit (q = 1), at the 
hypothetical tip of the wedge where the retardance is zero, the polarized light passes through 
the retarder and analyzer without hinderance. Moving along the slit, the retardance increases 
to the point where it becomes quarter-wave, and the light is converted to circularly polarized 
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light after the retarder, and half of the light transmits through the analyzer. As x, and the 
retardance, increase together, it reaches the point where there is half wave retardance. At 
that point the polarization is rotated 90° after the retarder, and none of the light transmits 
through the analyzer. At the same distance further along the slit, the retardance is full- 
wave, the light is completely transmitted, and the cycle is complete. Note that for typical 
birefringent materials, the spatial distance x corresponding to one wavelength of retardance 
will depend on the wavelength. In the absence of dispersion, the spatial modulation frequency 
is oc 1/A. 

Circularly polarized light {v = 1) is half transmitted at zero retardance (circularly polar- 
ized light passing through the analyzer). When the retardance reaches quarter- wave, the light 
becomes linearly polarized along the slit direction and all of the light passes the analyzer. 
When the retardance is half-wave, the sign of the circular polarization is flipped (v = — 1), 
and again, half of the light is transmitted through the analyzer. Hence, the modulation for 
v — 1 is similar to the modulation due to q = 1, but out of phase by a distance corresponding 
to one quarter wave of retardance. Light polarized linearly at 45° (it = 1), i.e., along the 
retarder fast axis, is unaffected by the variable retardance along the slit. However, if we pre- 
cede the birefringent wedge by an achromatic (or superachromatic) quarter wave retarder, 
with fast axis along the slit, then the circular polarization parameter V is interchanged with 
U [5]. Now V is unaffected by the variable retardance and causes no spatial modulation, 
while U causes spatial amplitude modulation, a quarter wave out of phase from Q. If the 
input Stokes vector is (J, Q, U, V), the output intensity I\\ in the spatial direction is 



The retardance <p maps onto the distance x along the slit according to <fi = 2ir(x/X), where 
the distance corresponding to a single wave of retardance change is X = A/(|n — n e \ tan£)}, £ 
is the wedge angle, and n Q and n e are the refractive indices for the o and e-beams, respectively. 
If the circular polarization V is desired rather than U, then the quarter wave retarder can 
be omitted. 

If a beam splitting polarizing prism (e.g. a Wollaston prism) is used as the polarization 
analyzer, then two versions of the spectra are obtained on the detector. The intensity I± of 
the orthogonally polarized beam is 



The difference of Eqs. (1) and (2), divided by their sum, assuming any transmission differ- 
ences have been removed, gives 



iH = 0.5(7 + Q cos 4> + U sin </>). 



(1) 



I ± = 0.5(7 - Q cos <j>-U sin 0). 



(2) 



V = - I±)/ (i|| + I±) = qcos(J) + usincj). 



(3) 
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Table 1 summarizes this and other algebraic expressions for the spatial modulation in sub- 
sequent configurations discussed below. An alternative way to express Eqs. (l)-(3) is 

i|| = 0.5/(1 +pcos(0- 2^)) 
I± = 0.5/(1 -pcos(0- 2^)) 

/' = pcos((j)-2ip), (4) 

where the position angle tp of linear polarization is given by ip = \taxr l {U/Q). From Eq. (4) 
it is apparent that the spatially modulated profile has an amplitude of modulation equal 
to the degree of polarization, and a (spatial) phase zero point that reveals the angle of 
polarization. 

It is important in using Eqs. (1) and (2) that there not be any signficant intensity variations 
along the slit on length scales of order X. However, in the "dual beam" version, Eq. (3), 
the total intensity / along the slit has been eliminated [5]. This potentially offers a means 
to retain some spatial resolution along the slit. For example, the image of a star can have 
very large intensity changes along a spectrograph slit, though its polarization is unchanged. 
Provided the extent of the image is sufficient to encode the sine and cosine terms in Eq. (3), 
we may derive its polarization even in the presence of quite strong intensity changes. In 
Eqs. (1) and (2), intrinsic intensity changes would be mixed with amplitude modulation 
produced by polarization. Hence, care needs to be taken in matching the spatial extent of 
the instrumental point spread function to the projected scale of the retardance variation. 
Use of Eq. (3) is more robust against this constraint. 

In practice, wedge components available off the shelf are relatively thick. This introduces 
a multiorder retarder effect, exploited in the SPEX concept [5]. A thick birefringent mate- 
rial, followed by analysis optics, such as those employed here, for a single location on the 
slit, yields spectral modulation, used to measure the polarization in [5]. Hence, using a sin- 
gle wedge, the resulting fringes from polarized light have a relatively pronounced "slope" 
because the retardance, implicitly, is varying as a function of both wavelength and spatial 
direction. From above, the complete expression for <fi is <p = 2irx(n — n e ) tan£/A. Therefore 
a constant phase <j), which defines the appearance of the fringes on the detector, occurs for 
x = [0/ (2ir(n e — n ) tan£)] A. The thickness of the wedge on the narrow edge effectively adds 
a constant to <p and hence increases the slope of x vs. A. 

In principle, analysis methods applied to the orthogonal dimension, discussed in Sec. 4, 
can remove this spectral modulation. The slope also translates to a constraint on the spectral 
resolution, as it must be sufficient to separate the fringes. 
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2.B. Double wedge 

An alternative strategy to alleviate the spectral resolution constraint, and at the same time 
simplify and render the analysis more robust, is to compound the first wedge with an iden- 
tical second wedge reversed in direction, with its fast axis orthogonal to the fast axis of the 
first wedge. The resulting optic, also used in Babinet-Soleil compensators, is convenient to 
work with (it is now a rectangular cuboid, see Fig. 3), and the reversed signs of birefrin- 
gence conspire to yield a zero retardance at the point where the two retardances are equal, 
expected to be close to the center of the optic. This is the configuration discussed by [2,7] 
in the context of imaging polarimetry. In imaging polarimetry, spectral modulation across 
the bandpass would be highly problematic, and this configuration goes some way towards 
eliminating it. [2] discusses other ways to mitigate birefringence dispersion. When applied 
to spectropolarimetry, as here, the effect is to remove the slope from the fringes and cause 
them to be approximately parallel in the spectral dimension, with spacing convergence to 
the blue, as a given amount of retardance corresponds to a greater number of wavelengths. 

Using Mueller matrix algebra, it can be shown that the amplitude modulated profiles for 
this configuration are, for the single beam, with a quarter wave plate used, as above, to 
impose a sensitivity to Stokes U rather than Stokes V: 

ill = 0.5(7 + Q cos 2<p + U sin 2<p) (5) 

and for the dual beam 

I' = q cos 2<p + u sin 2(f). (6) 

The essential characteristics of the compound wedge profiles are the same as the single wedge, 
though the spatial frequency is doubled and the multiorder retarder effect with wavelength is 
removed making the fringes essentially parallel to the dispersion direction. The double wedge 
also yields a quasi- zero order retarder, which has much smaller temperature dependence than 
a single wedge. 

Fig. 4 shows simulations of data frames obtained with 100% polarized light for each of the 
configurations discussed here, for Q, U and V polarized light separately. To illustrate the 
concepts of this section, Fig. 5 shows actual long-slit spectra obtained using one and two 
quartz birefringent wedges with a 3° wedge angle, installed in a simple slit spectrograph. 
We inserted the quartz wedges immediately after the entrance slit of the spectrograph, 
together with a quarter wave plate and analyzer as discussed. Fig. 5 shows fringes obtained 
when illuminated by linearly polarized light for both single wedge and double wedge. These 
correspond to configurations w and ww 1 defined in Table 1, though the wedges are not 
compounded in this test. The fringes are clearly visible when polarized light enters the 
spectrograph, and are not visible when unpolarized light is used. A more formal laboratory 
validation follows below in Sec. 6. 
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3. Data analysis Methods 

In equations (l)-(3), (5), and (6), the Stokes parameters are coefficients of orthogonal 
trigonometric functions. While Fourier methods could be used to retrieve these coefficients, 
we prefer a linear least squares solution. This enables us to generate formal error estimates, 
in addition to providing the Stokes coefficients in a straightforward fashion. With a least 
squares approach, one can also take the dispersion of birefringence into account in a more 
straightforward way than with Fourier methods. The general methods are described in the 
Appendix. If the spatial profile at a fixed wavelength is given by 

-^obs = Iic + Ql c + Uu c + Vv c (7) 

for a single beam, or 

/' = ~ I±)/(I\\ +I±) = q-q c + u-u c + v-v c (8) 

for the dual beam formalism, then the terms i c , q c , u c and v c are trigonometric functions 
whose coefficients are the Stokes parameters we seek, i.e., Eqs. (l)-(3), (5), (6), and others 
below. Since these trigonometric functions take on a set of known values at each Xj along 
the profile, standard methods can be applied to solve for their coefficients given an observed 
intensity profile I^^xi) or I'(xi), as described in the Appendix . The solution depends on 
the inverse of the curvature matrix, which is derived from products of the functions i c , q c , 
u c and v c . The solution is the inverse of the curvature matrix, multiplied by a vector derived 
from the observed profile and the same set of functions. The covariance matrix is also the 
inverse of the curvature matrix, and uncertainties on the Stokes parameters are taken as the 
square root of the diagonal terms of that matrix. 

Analytic expressions can be derived for the terms of the curvature matrix and its inverse. 
Replacing the summations by integrals over complete periods, it may be seen that in a 
formal mathematical sense the Stokes parameters are independent because the trigonometric 
functions are orthogonal for most of the methods described. That is, the integrals of their 
products become zero and the inverse of the curvature matrix is diagonal. Hence there is no 
formal covariance between the terms. In an instrumental realization this should translate to 
no formal cross-talk between the Stokes parameters. In the cases where there are off-diagonal 
terms in the curvature matrix inverse, we discuss this in the text. 

The description so far is idealized, however. In general, summations are not over complete 
periods, and departures of the wedge characteristics from those assumed, e.g. the exact value 
of the slope, and sampling errors can all yield off-diagonal terms. We discuss methods to deal 
with such issues in Sec. 5. These terms can generally be included in the analytic analysis, 
although the equations can become complicated depending on which tolerances are explored. 
An example can be found at the end of Sec. 4.B.3. 
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It will also be necessary to calibrate a given system. The goal of calibration can be either 
to determine the (hopefully) small correction factors to the analytic formulae, or else a 
completely empirical calibration approach can be adopted, as discussed in Sec. 5.B. 

To within an order of magnitude, we see in the analytic solutions, presented in the Ap- 
pendix, the expected \j\J N to t dependence for the sensitivity to which polarization can be 
measured, where N tot is the number of detected photons. 

4. Concept for Full Stokes Polarimetry 

The use of precision circular spectropolarimetry as a biosignature [3] requires accurate, sen- 
sitive measurement of the circular polarization spectrum. It is strongly preferred that all 
Stokes parameters be measured in order to better understand the physics involved and to 
guard against instrumental cross-talk, where one Stokes parameter is measured incorrectly, 
through instrumental effects, as another. Here, we describe two options to acquire full Stokes 
polarimetry data. The first, briefly described, is a brute force approach, where we simply 
place two versions of the configurations described above, next to one another. The second 
approach is to use an additional wedge or compound wedge with different fast axes and 
wedge angles to fully encode all Stokes parameters on a single data frame using a single 
optical bench. 

4- A. Two spectrograph slits 

In Sees. 2A and 2B, Eqs. (1) — (6), we showed how to encode two Stokes parameters simul- 
taneously. By dispensing with the quarter wave retarder in those configurations, e.g. using 
a single wedge single beam, we have 

7 Q = 0.5(1 + Q cos0 - V sin 0). 

If a second device is constructed with its wedge fast axis horizontal, 0°, and analyzer at 45° 
to the horizontal, then we have 

Iu = 0.5(7 + [/cos</> + Vsin0). 

These two configurations can be processed indendently, and typically \V\ <C \Q\, \U\. Fol- 
lowing the formalism of Eq. (4), the influence of a component of circular polarization v 
is to shift the phases of the spatially modulated profiles by amounts, in pixel space, of 
5x = X ta,n~ 1 (v/q)/(27i) for the first device and 5x = X tan~ 1 (f /m)/(27t) for the second de- 
vice, where X is the spatial period corresponding to one wavelength of retardance. For small 
v, the spatial shift of the profile is therefore Sx ps Xv/2nq and 5x ~ Xv/2nu, respecively. 
For example, if the source is 3% polarized, q = 0.03, if v = 10^ 3 (0.1% circularly polarized), 
and if X = 10 pixels, the spatial shift is 0.05 pixels. The accuracy to which the spatial shift 
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can be determined depends in principle on the total number of photons rather than the 
value of X, provided there are sufficient points within X to properly sample the profile (see 
Appendix). Systematic instrumental effects are also likely to enter; however, measurement 
to a precision of better than 0.01 pixels can be achieved in precision astrometry, and these 
accuracies ought to be feasible. 

4-B. Double wedges 

A concise solution to acquiring full Stokes spectropolarimetry may be obtained by using a 
more complex group of birefringent wedges. Our technique is analogous to the dual PEM 
polarimeters which encode the different Stokes parameters using different carrier frequencies. 
If we follow the first wedge or compound wedge by a second wedge that has twice the thickness 
gradient, and a fast axis oriented 45° to the first one, then the resultant intensity profiles 
encode the full Stokes information. The trigonometric functions involved are more complex 
than the simple sine and cosines of above, but they are relatively straightforward and still 
orthogonal. This configuration has the advantage of measuring the full Stokes parameters 
completely and simultaneously for a source, without minor viewing perspective differences 
in the double slit option (Sec. 4. A) and without needing to combine measurements from two 
essentially independent polarimeters. 

We describe three versions of this concept, though it is clear that a variety of permutations 
and options are available under the umbrella of this general approach. The first two versions, 
to illustrate the principles, are to place a single wedge with twice the thickness gradient par- 
allel and antiparallel, respectively, after a single wedge oriented as for the linear polarimetry 
application above. We conclude with a discussion of a double compound wedge, in which the 
second pair has twice the thickness gradient of the first. No quarter wave retarder is required. 
An example of a different permutation would be to reverse these two compound wedge pairs, 
which would result in a potentially more convenient choice of analyzer angle following the 
wedges. We defer discussion of such an option to a later paper and focus here on providing 
a proof of concept and demonstration of the approach. The naming convention used in the 
following subsection headings is described in the footnote to Table 1, which summarizes the 
terms in the equations for the amplitude modulation. Table 2 presents formulae for the uncer- 
tainties on the normalized Stokes parameters derived using the linear least squares methods 
of Sec. 3 and the Appendix. Table 3 gives the error estimates for the unnormalized Stokes 
parameters. Throughout, we assume standard (lcr), coverage factor k — 1, uncertainties. 

4.B.I. Two individual parallel wedges: wW 

In this version we assume the first wedge has a birefringence gradient, described by 4>(x), 
and fast axis at 45° to the horizontal, as above. We assume the second wedge has twice the 
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retardance 20(x) at location x and has fast axis aligned with the slit, 0°. Additionally, it is 
necessary to allow the analyzer angle to be set at angles other than zero degrees, in order 
not to lose dependence to the Stokes U parameter. We define the additional variable 9 as the 
angle of the transmission axis of the analyzer with respect to the horizontal or slit direction. 
In the dual beam configuration, the transmission axis of the second polarized beam is 9 + 90°. 
We retain the labelling In and I± for these two beams, respectively. In practice, there are 
likely to be constant offsets and a gradient of the second wedge not exactly twice that of 
the first wedge. These terms can generally be included in the analytic analysis of Sec. 3, but 
for clarity in introducing the concepts, we set them aside for now. It can be shown that the 
amplitude modulated profiles for this configuration are, for the single beam: 

ijl = 0.5(7 + Q(cos cos 29 + sin sin 20 sin 29) + U cos 20 sin 29 

+V(cos0sin20sin26) - sin cos 29)), (9) 

I± = 0.5(7 — Q(cos cos 29 + sin sin 2<j) sin 29) — U cos 20 sin 29 

-V(cos0sin20sin2# - sin cos 29)), (10) 

and hence for the dual beam, 

/' = g(cos cos 29 + sin sin 20 sin 29) + u cos 20 sin 29 

+v (cos sin 20 sin 29 — sin cos 29). (11) 

We see that if 9 = 0°, then the coefficient of U is zero and hence we cannot derive the 
value of U. There is no particular reason to have the analyzer at such an angle. For example, 
in Fig. 4 we use 9 = 45° for the double wedge configurations. However, the derived variance 
for each of the Stokes parameters is a function of 9. Depending on the application, it may 
be helpful to select 9 to minimize the variance of the estimate for Stokes V, since Stokes V 
is usually orders of magnitude smaller than the linear Stokes parameters. Below, Sec. 4.B.3 
we show an example. 

4.B.2. Two individual antiparallel wedges: wW 

If the second wedge is placed in the opposite direction to the first, so as to minimize the 
geometric angles of the two wedges together, we have phase gradients of 0(x) and (£ — 20(x)), 
respectively, for the two wedges, where ( is a constant, presumed unknown. It can be shown 
that the amplitude modulated profiles for this configuration are, for the single beam: 

J|! = 0.5(7 + Q(cos cos 29 + sin sin(C - 20) sin 29) + U cos(C - 20) sin 29 

+V(cos sin(C - 20) sin 29 - sin cos 29)), (12) 
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Jl = 0.5(7 - Q(cos cos 20 + sin sin(C - 20) sin 20) - U cos(C - 20) sin 29 

(cos sin(C - 20) sin 29 - sin cos 20)), (13) 



and for the dual beam 



I' = q cos cos 29 + sin sin(£ — 20) sin 29) — u cos(£ — 20) sin 29 

— v(cos0sin((^ — 20) sin26 l — sin cos 29). 



(14) 



In this case, where the second wedge is configured to run antiparallel to the first one, 
there is an off-diagonal term in the inverse curvature matrix B -1 (defined in the Appendix) 
containing the term (sin ( sin 49). (The same term appears if the wedges run parallel, but 
with a phase offset.) The inverse of the curvature matrix is the crucial mathematical entity 
in both solving for the Stokes parameters and estimating their variances. If the off-diagonal 
terms are zero then the trigonometric functions are orthogonal, and there are no formal 
dependencies of one Stokes parameter on the others or covariances between them. If there 
are off-diagonal terms, then there may be such dependencies and covariances. 

The off-diagonal term can be set to zero if the analyzer is placed at an angle 9 such that 
sin 49 is zero. Since the U coefficient also includes terms involving sin 29, which we do not 
want to be zero, the desired analyzer angle to eliminate cross- dependencies is 9 = 45°. There 
may be a trade choice based on the optimization required, between making the off-diagonal 
terms of the inverse curvature matrix zero and minimizing the variance of a particular Stokes 
parameter, which can occur at a different analyzer angle. 

4.B.3. Two compound wedge pairs: ww'WW' 

By analogy with the compound double wedge in the linear polarimetry example, described 
in Sec. 2.B, it is possible to use two compound double wedges to provide zero retardance 
in the center of the optic. The first wedge pair has one gradient, as for the optic used in 
Sec. 2.B above (W in the notation of Table 1). Each element of the second wedge pair has 
double the thickness gradient. The first has its fast axis along the slit, while the second has 
its fast axis orthogonal to it, so that the retardances are £ — 2<p(x) and £ + 20(x), respectively. 
Again, the compound device with four wedges is rectangular in shape. It can be shown that 
the amplitude modulated profiles for this configuration are, for the single beam, 



J|l = 0.5(7 + Q(cos20cos20 + sin 20 sin 40 sin 20) + U cos 40 sin 20 

+V"(cos 20 sin 20 sin 20 - sin 20 cos 20)), 



(15) 



I± = 0.5(7 - Q(cos20cos20 + sin 20 sin 40 sin 20) - £/ C os40sin20 

— V (cos 20 sin 20 sin 20 — sin 20 cos 20)), 



(16) 
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and for the dual beam 

/' = g(cos 20 cos 29 + sin 20 sin 40 sin 29) — u cos 40 sin 29 

— (cos 20 sin 40 sin 2$ — sin20cos2#). (17) 

These functions are orthogonal to one another. Fig. 6 shows the derived variance as a 
function of the analyzer angle, with analytic formulae presented in Table 2. The minimum 
variance for Stokes V, which typically has the smallest value of the Stokes parameters, is 
given at analyzer angle with tan 4$ = —2, which implies 9 ~ 74.1°. The minimum variance 
for Stokes Q is at tan 49 = +2, i.e., 9 ~ 15.9°, and that for U is at 9 — 45°. The value of 
the minimum variance is, for q and v, a(q,v) ~ 1.24/ y/N to t, which represents the photon 
counting limit for such a device and which is quite close to the canonical 1/ a/ N to t value. 

If the centers of the two compound wedges are misaligned by a spatial distance s, then 
there is an additional term in the above equations, which we characterize by a phase offset in 
the second wedge pair. That is, their retardances are given by £ — 20(a;) — a and £ + 20(x) + a, 
respectivel, where a = 4-ks/X. In this case, the amplitude modulation expressions are 

J|l = 0.5(7 + Q(cos 20 cos 29 + sin 20 sin(2a + 40) sin 29) + U cos(2a + 40) sin 29 

+K(cos 20 sin(2a + 40) sin 29 - sin 20 cos 29)), (18) 

I ± = 0.5(7 - Q(cos 20 cos 29 + sin 20 sin(2a + 40) sin 29) - U cos(2a + 40) sin 29 

-V{cos 20 sin(2a + 40) sin 29 - sin 20 cos 29)), (19) 

and for the dual beam 

/' = g(cos 20 cos 29 + sin 20 sin(2a + 40) sin 29) - u cos(2a + 40) sin 29 

-v(cos 20 sin(2a + 40) sin 29 - sin 20 cos 29) . (20) 

Following through the least squares analysis, it can be seen that the miscentering intro- 
duces an off-diagonal term in the inverse of the curvature matrix, correlating the errors in 
Q and V, and again involving sin 49. Hence, this off-diagonal term can be set to zero by 
setting 9 = 45°. If this is not done, or the tolerances do not allow it, then the presence of a 
covariance term in and of itself does not bias the result. The random errors on Q correlate 
with the errors on V. The actual values of Q and V do not correlate, though, provided that 
the covariance term is known. In principle the covariance term does contribute to the overall 
error budget, although assuming independence yields a good description of the actual ob- 
served variance in Monte Carlo simulations, shown in Fig. 7. In the case where the term is 
present, but known, we still correctly derive the Stokes parameters without cross-talk. 
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Cross-talk occurs if the covariance is not characterized correctly. For the example here of 
two miscentered wedge pairs, if we characterize the miscentering by a/X, the ratio of the 
offcenter distance to the spatial distance corresponding to one wavelength in the first wedge 
pair, then the true value of Stokes v, say, is v true = sy^B" 1 ]^, + sy^B' 1 ^, following the 
notation of the Appendix. If we incorrectly derive v es t = sy4[B~ l ] vv , ignoring the cross-term, 
then v true — v est = s?/i[.B _1 ] g „ = g[_B _1 ] g „/[i? _1 ]gq. A formal tolerance analysis can be carried 
out, and to completely ignore this term, while restricting the cross talk < 10~ 3 , requires 9 to 
be within < 0.1° of 45° if the miscentering is such as to maximize the cross talk. If needed, 
there are two options to relax this constraint within the context of an analytical approach: 
(i) do not ignore the term! and (ii) increase the spatial scale X to ease the requirement on 
s/X. Empirical calibration approaches are also possible as discussed below. 

5. Practical Implementation, Requirements, Sensitivities and Calibration 

Wedge retarders, such as those discussed here may be custom manufactured. However, 
testable quality versions are available off-the-shelf under the guise of depolarizers or scram- 
blers. Quartz or calcite provide plausible birefringent materials. The ordinary and extraor- 
dinary refractive indices and birefringences are, for quartz, n Q = 1.5384, n e = 1.5473, 
n e — n = 0.0089, respectively, and for calcite n Q = 1.647, n e = 1.480, n e — n = —0.167. 
For example, if the retarders have a wedge gradient of 3°, then at 500 nm the retardance 
increases by one wavelength over a distance X of ~ 1.1 mm for quartz and 0.06 mm for 
calcite. If these scales are projected 1:1 onto a detector with 5 /jm pixels, then these values 
correspond to X m 214 and 11 pixels for quartz and calcite, respectively. 

5. A. An empirical calibration approach 

In a similar fashion to the analysis described in Sec. 4.B.3, the least squares methods lend 
themselves to formal tolerance analyses, as well as parameter estimation. However, a compre- 
hensive analysis of all plausibly relevant parameters is impractical and premature. Instead, 
we consider an alternative approach, which is purely empirical. For a very general set of 
optical component characteristics, the generic versions of the amplitude modulated intensity 
profile, Eqs. (7) and (8), are valid, even if the exact forms of the functions i c , q c , u c and 
v c are not known. If we present the system in turn with unpolarized light, and then 100% 
polarized light oriented in the Q direction, the U direction, and finally, with 100% circularly 
polarized light, then the empirical response gives the functions i c , q c , u c and v c . These em- 
pirically derived functions can then be used to numerically derive the curvature matrix and 
its inverse. As in Sec. 4.B.3, the presence of covariance terms in the matrices by itself does 
not invalidate the approach, because, in principle, with a high quality set of calibration ob- 
servations (which would appear much like the examples shown in Fig. 4), they are implicitly 
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known. It is only when the terms are not correctly accounted for, that problems may arise; 
that is, if the calibration sources are not of high enough quality. It is also likely that the 
variance will be a function of analyzer angle, as above. Hence empirical versions of Fig. 7 
may be useful to explore trade space. 

We expect, in general, the empirical calibration method to be a very important approach, 
potentially offering the best strategy for deriving the Stokes parameters accurately. However 
it remains to be seen to what degree the required tolerances, or calibration stability, limit 
the performance of these devices in practice. Our intent is to develop additional laboratory 
experience to understand these issues. 

5.B. Sensitivities 

Component birefringence will depend on temperature. This can be mitigated by (i) stabi- 
lizing the temperature, (ii) continuous calibration, and (iii) compounding carefully chosen 
materials. To obtain an idea of the order of magnitude of the temperature sensitivity, we 
use the temperature dependent formula for the birefringence B = n e — n of quartz, given 
by [21], which yields (dB/dT)/B « 1.2 x 1(T 4 for the range T = °C to 25 °C. Since the 
wedge should keep its shape under expansion or contraction, only the birefringence term 
matters and imprints itself as an identical fractional change e = 1.2 x 10~ 4 dT on the spatial 
wavelength X. It can be shown that the resulting fractional change in q for the dual beam 
single wedge example is one half the fractional change in B, for an input beam consisting 
of purely q, and the spurious cross-talk v! into the other Stokes parameter u is u' = 7req. 
Hence a 1 °C temperature change would result in a spurious polarization (cross talk) of 
v! ~ 3 x 10~ 5 for a 10% linearly polarized source. 

The presence of two refractive indices in a wedge-shaped optic will cause a prismatic 
separation of the orthogonally polarized beams. The magnitude of this effect will depend 
on the details of the optical system designed. In our laboratory testing, this issue was not 
significant. 

If the beam incident on a single wedge has significant convergence, then the retardance 
seen by light at different angles differs by approximately 1 / cos ( where ( is the angle to the 
normal. For this retardance difference to be less than A/10, say, the /-number of the beam 
needs to be slower than / /9 for a 1 mm thickness quartz and //40 for 1 mm thickness calcite. 
The convergence requirement scales as the inverse square root of the thickness. 

We carried out a number of Monte-Carlo simulations and found that the retrieved solutions 
are close to the theoretical expressions given in the tables provided that (i) at least one full 
period is sampled and (ii) there are sufficent sampling points within the full cycle to properly 
sample the highest frequency component. 

In cases where the fringes have significant slope, as described above, it will be necessary 
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for the spectral resolution to be sufficient to separate the fringes. If we require the retardance 
change A<j)/(2ir) < 0.1, this is (d<f)/d\)A\ < 0.1, where <fi = 2irL(n e — n )/\, L the thickness. 
Hence the required spectral resolution is R = -A- = 10L(n e — n )/A. For 2 mm thick quartz 
this is R > 360 and for 1 mm calcite, R > 3330. For the compounded wedge optics which 
straighten the fringes, this constraint is relaxed to the point of being essentially irrelevant. 

6. Laboratory Validation 

We established an optical test bench to allow us to provide an empirical proof of concept 
for the approach presented in this paper, and to demonstrate that the device functions as a 
polarimeter. The optical test bench is shown and illustrated in Fig. 8. 

Light sources, either halogen continuum white light, or a variety of line lamps for wave- 
length calibration, illuminated an integrating sphere's entrance port. The light emerging 
from a separate port at right angles to the first is expected to be unpolarized. However, to 
ensure an unpolarized source and to provide a uniform location for our measurements, the 
emergent light was directed to fall onto an opal diffusing screen. Following the opal screen, 
and located close to the screen, the light optionally encountered polarizing elements (Q, U, 
or V) for calibration, samples to measure the polarized transmission spectrum, or nothing, 
to provide an intensity reference spectrum. 

The spectrograph consisted of a slit, a collimator, a transmission grating, and a camera. 
The slit was 125 /im by 1 cm and was located at one focus of an //1.4 50 mm Nikon 
collimating lens. 1 The transmission diffraction grating was ThorLabs Part # GT50-06V with 
600 grooves/mm. A second //1.4 50 mm Nikon lens imaged the spectrum onto a Quantum 
Scientific Imaging 683 8 Mpix cooled CCD Camera with 3326 x 2504 pixels of size 5.4 micron. 

For the polarization analyzer, we placed a Meadowlark precision linear polarizer, 
DPM100VIS between the spectrograph slit and the collimating lens, using a rotary mount 
to allow adjustment of the analyzer angle. The extinction ratio across the visible spectrum 
for this polarizer is ~ 10 4 , and exceeds 100:1 from approximately 375 nm to 725 nm. The 
birefringent wedges were placed on a platform next to the spectrograph slit between the slit 
and the source (i.e. before the light enters the spectrograph). The optional quarter wave 
retarder was mounted, when used, between these wedges and the source and was an achro- 
matic Meadowlark quarter wave retarder AQM- 100-545, which is effective between 450 nm 
and 630 nm. The birefringent wedges themselves were a mix of customized and off-the-shelf 
quartz scramblers from Karl-Lambrecht Corporation. Off-the-shelf wedges had a pitch angle 

1 Certain commercial equipment, instruments, or materials are identified in this paper in order to specify 
the experimental procedure adequately. Such identification is not intended to imply recommendation or 
endorsement by the Space Telescope Science Institute, the National Institute of Standards and Technology, 
or the Universiteit Leiden nor is it intended to imply that the materials or equipment identified are necessarily 
the best available for the purpose. 
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of pa 3°, and the customized pieces used a pitch angle of ~ 6°. 

The distance between the source and the spectrograph slit was approximately 0.6 m, and 
the entire system was contained within a series of light-tight baffles and boxes to eliminate 
stray light. Locating the birefringent wedges externally to the spectrograph not only allows 
us to ignore possible polarization due to the slit, but also allowed easy access for switching 
configurations, and permitted us to focus the spectrograph onto a single slit position. 

The first exercise was to attempt to reproduce the appearance of the theoretical data 
frames for the various configurations presented in Fig. 4. Fig. 9 shows the results. To do this, 
we used the white halogen continuum source in combination with linearly and circularly 
polarizing filters. For the first four configurations (first four columns in Fig. 9), we used a set 
of 60 mm astronomical polarizing filters that utilize HN38 Polaroid mounted in a magnesium 
fluoride substrate to approximate 100% linearly polarized light (first two rows of the first 
four columns in Fig. 9). We used a cholesteric liquid crystal technology (CLC) filter to 
approximate 100% circularly polarized light (third row of the first four columns in Fig. 9). 
For the final configuration (fifth column of Fig. 9), we used a polarization state generator 
that utilized a precision linear polarizer in combination with a Fresnel rhomb [22] that is 
capable of producing close to 100% polarized light anywhere on the Poincare sphere. We also 
obtained data frames without any polarizing optics to provide a "flat-field" reference. The 
spectral scale was wavelength calibrated using argon and mercury line lamps. We defined 
the slit direction to correspond to Q and 45° to the slit to correspond to U. 

By comparing Figs. 4 and 9, it is apparent that the empirical data reproduce the qualitative 
expectations extremely well. Differences in detail can be attributed to different absolute 
thicknesses of the wedges (in practice the 3° and 6° wedges had very different thicknesses) 
relative to one another and to the theoretical model, and to essentially random centering of 
the crossover points for wedge pairs relative to the slit and to one another. Care was taken 
with the sign convention and parity of the wedges to reproduce the directions for fast axes 
and wedge gradients used in the models. The empirical wavelength range shown corresponds 
to 550 nm to 700 nm, and the slit height to 2 mm in Fig. 9. We consider the data obtained 
in this exercise to be fully consistent with the theoretical expectations, given the practical 
uncertainties described. 

The second exercise was to test the linear polarization mode described in Sec. 2B. We 
placed a compound 3° wedge pair with the fast axes crossed, running corner to corner at 
±45° to the spectrograph slit, in front of the spectrograph slit. An achromatic quarter wave 
retarder was placed upstream of the wedge pair, with its fast axis oriented 0° to the slit. 
Calibration measurements were taken using the standard suite of polarizing filters, and the 
resulting data frames were used to derive empirical "coefficient" frames for use in the least- 
squares retrieval procedure of Sec. 3 and Sec. 5A, with additional theoretical analysis in 
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the Appendix. We placed an Schott BG18 colored glass filter (broad bandpass with peak 
transmission near ~ 510 nm) close to the source and orthogonal to the beam and took 
a single data frame. We then rotated the filter about the vertical axis by approximately 
30° and took a second data frame. Using these single data frames in conjunction with the 
procedures described in Sec. 3 and Sec. 5A, we derived polarization and intensity spectra, 
shown in Fig. 10. The average linear polarization for the region 500 nm to 550 nm was 0.39% 
and 6.2% with standard deviations 0.02% and 0.35% respectively. For a uniform glass sheet 
of refractive index n = 1.5 in air, we expect transmitted light to be polarized 6.3% for an 
incident angle of 30°. We emphasize that these results were obtained using a single data 
frame with no moving parts once the calibration data had been acquired. 

The third exercise was to demonstrate our ability to obtain full Stokes polarimetry from 
a single data frame. For this, we used the compound 3° wedge pair used in the previous 
exercise, followed by a compound 6° wedge pair with fast axes at 0° and 90°. We removed 
the quarter wave retarder. As an interesting source of circularly polarized light, we used a 
pair of plastic 3D cinema glasses. These glasses comprise a quarter wave sheet and polarizing 
sheet in combination. Used as viewers, they either transmit or extinguish circularly polarized 
light, depending on its sign. But, in reverse, they produce circularly polarized light of opposite 
signs for each "eye." We took single data frames through each eye in turn and processed 
the data frames according to the methodology of Sec. 3 and Sec. 5A, using the empirical 
suite of calibration data frames. The results are shown in Fig. 11. The method produced 
excellent results on these sources, yielding the expected extremely high level of polarization, 
with opposite sign for the two eyes. We also show the derived linear polarization, to show 
that we are also measuring full the Stokes vector in these single data frames. 

We defer attempts to carry out precision polarimetry given the rudimentary nature of our 
optical bench coupled to imperfect calibration source availability. However, we believe that 
these preliminary results are extremely encouraging and satisfy our desire to provide a proof 
of concept in the laboratory and to demonstrate that the desired polarimetric information 
can be retrieved from single data frames. 

7. Conclusions 

We have described an approach to polarization measurement that uses no moving parts and 
that relies on simple, robust optical components. Either linear polarimetry or full Stokes 
polarimetry can be carried out. The method depends on the use of an area detector, such 
as a CCD, with the light spread across a region of the detector. If the system can be made 
photon-limited in sensitivity, this spreading of the light improves the polarimetry, since 
typical CCD well-depths are only of order 10 5 . With modest spreading of the light, a single 
photon limited frame should be able to reach precision of order 10~ 4 in polarization. The 
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influence of departures from ideal circumstances still remains largely to be explored. Hence, 
we do not know at this stage whether this approach will be able to achieve extremely high 
accuracy. However the robustness and simplicity of the components involved offers cause for 
optimism. Other approaches, such as the spectral modulation method for linear polarimetry 
[5], offer alternative methods for static polarimetry in hostile environments. Compared to 
that approach, the methods presented here yield a cleaner separation of the spectroscopy and 
polarimetry, at the expense of additional detector surface area requirements. The methods 
may be applied in the UV or IR as well as in the visible wavelength range. 

Since the entire polarization information is contained within a single data frame, the 
method is well-suited to measuring the polarization of transient sources and scenes where the 
polarimeter and target are in rapid relative motion. Since the optics are robust, simple and 
require no moving parts, we anticipate that these methods will prove useful for application 
in space. 
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Appendix: General Linear Least Squares Methods 

We follow Bevington [23, 24] and let the general problem to be solved be 

y(xi) = ai c (xi) + bq c (xi) + cu c (xi) + dv c (xi), 

where measurements jji, either the intensity / in the single beam case or (Jy — /_!_)/( Jy + 1±) 
in the dual beam case, are made at points Xi and j/j = y(xi) + 6j, with y the true underlying 
value and e,; its error (assumed random, independent) at location x^. The terms i c , q c , u c , 
and v c are trigonometric functions that encode the Stokes parameters /, Q, U, and V or q, 
u, and v, and their coefficients a, b, c, and d are the Stokes parameters to be derived. The 
mapping and specific functions depend on the chosen configuration, but all configurations 
discussed here can be expressed in this way. Sometimes the functions are identically zero, 
implying no sensitivity to that parameter. 
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The x 2 function is then 

N N N 

X 2 = ^2 e2 i/ a i =^2—[yi- y( x i)? = ^2—[Vi- aic(xi) - bq c (xi) ~ cu c (xi) - dv c (x t )} 2 , 

i=i i=i 0i i=i 0i 

and to solve, we set the partial derivatives of x 2 with respect to each of a, b, c, and d equal 
to zero: 

d X 2 



da 
dx 2 

db 

dx 2 
dc 

dx 2 
3d 



= -2 E \% c (yi - ai c - bq c - cu c - dv c ) , 

i 

= -2 E \q c (yi - ai c - bq c - cu c - dv c ) , 

i 

= -2 E \u c (yi - ai c - bq c - cu c - dv c ) 

i 

: = -2 j?Vc (Vi ~ aic ~ bq c - cu c - dv c ) 



We require the curvature matrix B and summation vector s 



.</• 



B = 



( E ^i 2 E ~^icq c E ~^ic u c E ^icVc \ 

i i i i 

Y^cUc YjtQcUc Yh U l Y,7? U cVc 

i i i i 

V E jticVc E jlVcVc E ^ U cVc E I 



respectively. With this terminology, the least squares equations become 



Sy = B- 



b 

c 

\dj 



We solve for the vector (a, b, c, d), 



a = (a, b, c, d) — B 1 ■ s y . 

Following standard procedures, e.g. [23], [24], ignoring covariances, the uncertainties on these 
parameters are 

= B ii 

where dj represents a, b, c, or d, and [-B _1 ]n is the corresponding diagonal term of B _1 . 
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Our application is precision polarimetry, for which it is presumed (i) the degree of polar- 
ization is small, and (ii) light levels are relatively high. Hence, the intensity across the spatial 
segment is approximately constant and obeys Poisson counting statistics. That is, we assume 
the uncertainty is the same for each bin, <Tj = a = (Ny/nx)^ 1 ^ 2 ^ where JV 7 is the total number 
of detected photons, and nx is the number of bins across which the photons are distributed, 
i.e., the number of x sampling points. For cases where the trigonometric functions that are 
embodied by i c , q c , u c , v c are orthogonal (we approximate the summations by integrals over 
integer numbers of periods), B is diagonal. Hence its inverse is also diagonal, and the Stokes 
parameter solutions are independent of one another. 

We choose as a simple example the single wedge with its fast axis at 45° to the slit direction, 
which, in turn, defines the direction for Stokes Q. In general, we use Mueller matrix algebra 
to solve for the system. As in Sec. 2. A above, Eq. (1) gives the expression for the intensity at 
points xf. yi = I(xi) = 0.5(7 — Q cos <f>i — U sin sin 0j) where 0; = 2n(xi/X). In the formalism 
above, y { = ai c (xi) + bq c (xi) + cu c (xi) + dv c (xi) so i c (xi) = 0.5, g c (x;) = 0.5 cos(27rXi/X), 
u c {xi) = 0.5 sm(27TXi/X) and v c (xi) = 0. In the absence of V the matrix B is reduced to the 
3x3 matrix 



B 



a' 



( X* c 2 

E»c?c 



E?c«c 



The summations run across nx pixels. Applying the expressions for the g c , u c , and v c , we 
have 



B 



4cr 5 



/ nx J2 cos (^ nx i/-^) E sm (27T^i/A) 

X cos(27rxi/A) X cos2 ( 2 W^O J2 cos ( 27r Xi/X)sm(2TTXi/X) 

\ sm(2TVXi/X) cos(2vrx i /A) sin(2vrx i /X) £ sin 2 (27rx i /A) 



/ 



We assume the summations cover an integer number of periods and approximate the sums 
using integrals, using f(x)dx ~ Ax Y^T f{ x i)i where Ax = X/nx (the width of a bin in 
x). Hence, Ylfi x i) ~ x / f{ x )d%- Thus, it can be shown for this example that 

/ 2 

nx 

B = — 1 



and 



B 



8(T 2 



4a 2 



\ 



o\ 



1 / 



nx 



f 1 





\ 


2/ 



Now, we can go back to the expression for s y and solve for the Stokes parame- 



ters, noting that X Vi 



X 



tot 1 



the total number of photons collected, and s t 
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2^2 E V ; f Yl C0S (2vrXj/X)?/j, sin(27rx i /X)?/j], omitting the zero V term. Hence the solu- 
tions from (a, 6, c, d) = B -1 ■ s y are 

I = a = (ff)(^) =^* = 2(y), 

\ 2a z J \ nx J nx >«* ' ' 

Q = & = (^)E2/iCos(27r^/X), 
^ = c = (^E^sin^/X). 

Similarly, we can derive the uncertainties of the Stokes parameters. The uncertainty o is 
given by a 2 = N tot /nx and is assumed to be constant. Hence, reading directly from the 
expression for B _1 , we have 

a(U)=a c 



nx 



The uncertainties in the normalized Stokes parameters q and u are 

a[q) = ^ = si- 

Ignoring bias terms, it follows that the uncertainty on the degree of polarization is 



cr(p) 



N, 



lot 



The expressions for the trigonometric functions i c , q c , u c , v c depend on the configuration 
and whether a dual beam formalism is adopted or not. We derived the expressions for these 
functions using Mueller matrix algebra for a selection of configurations, as presented in 
Table 1. In a similar fashion to this example, though with more complex manipulations, we 
can analytically invert the corresponding curvature matrix to derive both the solution and 
the uncertainty estimates, taking only the diagonal terms as the uncertainty. There are cases 
where the off-diagonal terms of B are non-zero, as discussed in the text. 
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Table 1. Coefficients of Stokes parameters for different wedge configurations 



Wcdgcs a 


Beam 






Qc 


u c 






qw 


single 


0.5 




0.5 cos 


0.5 sin 








dual 






COS 


sin 






qww' 


single 


0.5 




0.5 cos 20 


0.5 sin20 








dual 






cos 20 


sin 20 






wW 


single 


0.5 


0.5(cos 0cos 26 H 


- sin sin 20 sin 26) 


O.5cos20sin20 


0.5(cos sin 20 sin 26 - 


- sin cos 26) 




dual 




cos cos 26 


+ sin sin 20 sin 26 


cos 20 sin 26 


cos sin 20 sin 26 


— sin cos 26 


wW' 


single 


0.5 


0.5(cos cos 20 + sin q 


>sin(C - 20) sin 26) 


0.5 cos(C - 20) sin 20 


0.5(cos 0sin(C — 20) sin 26 - 


- sin cos 26) 




dual 




cos cos 26 + sin 


sin(C — 20) sin 26 


cos{C - 20) sin 26 


cos sin{£ — 20) sin 26 


— sin cos 26 


ww WW 1 


single 


0.5 


0.5(cos 20 cos 26 + 


sin 20 sin 40 sin 26) 


0.5 cos 40 sin 26 


0.5(cos 20 sin 40 sin 26 — 


sin 20 cos 26) 




dual 




cos 20 cos 26 - 


- sin 20 sin 40 sin 26 


cos40 sin 26 


cos 20 sin 40 sin 26 - 


- sin 20 cos 26 



"Notation for wedge configurations: q denotes optional quarter wave retarder, w denotes thickness gradient 
4> w = lixxjX and W denotes twice the thickness gradient <pw = Airx/X. A primed symbol denotes antipar- 
allel wedge direction relative to unprimed. Wedges w have fast axis at 45° to the slit, w' at —45°; wedges W 
have fast axis at 0°, except in combination WW, when they are at Oo and 90° respectively. 



Table 2. Error estimates for normalized Stokes parameters for different wedge 
configurations 



Wedges 


Beam 


qw 


single 




dual 


qww' 


single 




dual 


wW 


single 




dual 


wW' 


single 




dual 


ww WW' 


single 




dual 



(2/iV tot ) 1/2 
(2/N to t) 1/2 
(2/Ntot) 1 ' 2 
(2/Ntot) 1 ' 2 
2(2 /N tot) 1 ' ' 2 (3 + cos 49 + 2 sin 49) 
2(2 /N tot) 1 ' " 2 (3 + cos 49 + 2 sin 49) 
4/Jv l/2 / 3 + cos 49 + 2 cos C sin 49 ^ 
^' lv tot ^ 15+12 cos 49 + 5 cos 88 ) 
4/Jv l/2 / 3 + cos 49 + 2 cos C sin 4» W 2 
*'"tot \ 15+12 cos 19 + 5 cos 89 J 
2(2 /N tot) 1 /2 (3 + cos 49 + 2 sin 49) ~ 1/2 
2(2 /N tot) 1 ' 2 (3 + cos 49 + 2sin49)~ 1 '' 2 



o-(u) 



(2/Ntot) 1 ' 2 
(2/Ntot) 1 ' 2 
{2/Ntot) 1 ' 2 
(2/Ntot) 1 ' 2 



cr(v) 



1/2 
1/2 
1/2 



(2/Ntot) 1 ' 2 / 


sin 29 


(2/Ntot) 1 ' 2 / 


sin 29 


(2/Ntot) 1 ' 2 / 


sin 29 


(2/Ntot) 1 ' 2 / 


sin 29 


(2/Ntot) 1 ' 2 / 


sin 29 


(2/Ntot) 1 ' 2 / 


sin 29 



2(2/JV tot ) 1 / 2 (3 + cos 49 - 2sin49)~ 1 / 2 
2(2/JVtot) 1/2 (3 + cos 49 - 2 sin 49) ~ 1/2 
S49-2 co s C s 



4/N 1 ' 2 (3+c 

*/"tot \ 15 + 12 cos 49 + 5 c 



1/2 / 3 + cos 49 — 2 cos ( 
tot \ 15 + 12 cos 49 + 5 



^nV 2 

89 ) 

,49\l/2 



( 

2(2/JViot) 1 / 2 (3 + cos 40 — 2 sin 46) 1 / 2 
2(2/ Ntot) 1 / ' 2 (3 + cos 49 - 2 sin 49) _1 / 2 



Table 3. Error estimates for unnormalized Stokes parameters for different 
wedge configurations 



Wcgcs 


Beam 


o-(I) 


qw 


single 
dual 




qww' 


single 
dual 




wW 


single 
dual 




wW' 


single 
dual 




ww'WW' 


single 
dual 





"(Q) 



a(U) 



a(V) 



2(2Ntot) 1/2 /nx 
2(2N t ot) 1/2 /nx 



(Ntot S -± 



,48 V/ 2 



+ 12 cos 49 + 5 c 



2(2N t ot) 1/2 /nx 
2(2N t ot) 1/2 /n* 
2(2N t ot) 1/2 /(nx\ sin29|) 
2(2JV tot ) 1 / 2 /( na: | s i n 2e|) 



[Ntot 



3+cos 49-2 cos C sin 49 ^ V 2 
15 + 12 cos 49 + 5 cos 88 J 
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Fig. 1. Illustration of the optical bench layout. Light enters a spatially elon- 
gated slit, passes through a birefringent wedge or wedges and polarization an- 
alyzer, and then enters a conventional long-slit spectrograph. A quarter wave 
retarder (not shown) may be inserted before the wedge. 
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Fig. 2. Concept — the polarization optics imprint an amplitude modulation 
on the dimension orthogonal to the dispersion direction of the spectrograph. 
Typically, this direction corresponds to the spatial dimension along the slit. 
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Fig. 3. Illustration of the compound birefringent wedges. A single compound 
optic may be used for linear polarimetry (left) and double (both) for full Stokes. 
The fast axes run at ±45° in the left optic, horizontal and vertical in the right 
optic, and the slit direction is horizontal or vertical. 
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Fig. 4. Example theoretical data frames for each of the configurations dis- 
cussed in the text when viewed with 100% polarized light. Each panel has x 
running horizontally and wavelength vertically, increasing up. Parameters cor- 
respond to 2 mm in x of a 3° quartz wedge set, running from 450 nm to 750 nm. 
Top row shows 100% Stokes Q, middle row 100% Stokes U and bottom row 
100% Stokes V. Left to right, in the notation of Tables 1-3, the configurations 
are qw, qww', wW, wW, and ww'WW. Note that if the quarter wave retarder 
were omitted in the first two columns, U and V, which show no sensitivity with 
the quarter wave retarder, would be interchanged. For the first two columns, 
the fast axis is set at 0° and for the remaining three at 45° (see text). 
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Fig. 5. Actual spectra obtained with a preliminary test optical bench: (upper) 
a spectrum obtained with unpolarized light, (center) a spectrum with one 
quartz birefringent wedge and 100% linearly polarized light, and (lower) a 
spectrum obtained with two quartz birefringent wedges reversed as in the 
manner of the compound optics with 100% linearly polarized light. 
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Fig. 6. The uncertainties on normalized Stokes parameters for double wedge 
pair configuration, Sec. 4.B.3. Blue is Stokes q, green is u, and red is v. The 
vertical lines indicate the positions of the minima for q and v, the horizontal 
lines indicates 1 / y/ N tot (dotted), and the analytically-determined minimum 
1.24 x higher (dashed). 
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Fig. 7. The uncertainties on normalized Stokes parameters for double wedge 
pair configuration, Sec. 4.B.3 with miscentering of order l/32x the spatial 
distance of one wavelength of retardance. Blue is Stokes q, green is u, and 
red is v. The vertical lines indicate the positions of the analyzer angles that 
have no formal covariance, which is independent of the miscentering. Smooth 
lines through the simulated data are the analytic solutions ignoring covariance 
terms, while the plus signs are the results of Monte-Carlo simulations. The 
horizontal lines are as in the previous figure. 
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Fig. 8. Laboratory optical bench layout as implemented. Light enters an in- 
tegrating sphere, and illuminates a diffusing screen on exit. It then passes 
through calibration or sample components, before entering the long slit spec- 
trograph with its associated polarization components as described in the text. 
The upper panel shows the actual optical bench with baffles and boxes re- 
moved for visibility. The integrating sphere is to the left, and the camera to 
the right. 
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Fig. 9. Example data frames for each of the configurations discussed in the 
text when viewed with m 100% polarized light, obtained in the laboratory. 
The rows and columns correspond to those shown in Fig. 4. 
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Fig. 10. Retrieved polarization curves for a BG18 colored glass filter presented 
orthogonal to the beam, blue, and at an angle tilted by ~ 30° to orthogonal, 
green, observed using the wW configuration. At right angles, we expect no 
polarization, and inclined at 30°, approximately 6.3%, consistent with the least 
squares retrieval. The black curves show arbitrarily normalized throughputs 
for the two configurations (solid, orthogonal and dotted, inclined) derived from 
the data, serving to illustrate that we also obtain full Stokes I spectroscopy 
using these methods. 
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Fig. 11. Retrieved circular polarization for a pair of polarizing cinema 3D 
glasses, expected to exhibit 100% Stokes V left and right circularly polarized 
light for the left and right eyes, measured using the ww'WW configuration. 
The retrieval is consistent with expectations. For completeness, and to illus- 
trate that we obtain full Stokes polarimetry from a single data frame, the 
dashed lines show the retrieved degree of linear polarization. 
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